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Abstract 

The Monte Carlo generator EKHARA allows to simulate the process e + e~ — > e + e~n° at energies of meson factories. 
We review the motivation, describe the current status of the generator, and discuss features, which will appear in the 
forthcoming releases: the two-photon production of 77 and v[ mesons. 
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1. Introduction 

The two-photon transition form factors of the light 
pseudoscalar mesons (P - 7r°, 77, 77') are being inten- 
sively discussed during the last years. Recently the 
BaBar experiment has provided us with new informa- 
tion in the region where the virtuality of the first photon 
is high (g 2 ) and that of the second one is small (and 
both photons are spacelike) EE). The process, which 
gives an access to the given kinematic region, is e + e~ — > 
e + e~ f P in the "single-tag" mode, see Fig. [T] (left). Within 
the range of the previous experiments (CELLO 0, 
CLEO El) there was a fair agreement with the theory. 
The results of BaBar on the pion transition form fac- 
tor O demonstrated an unexpected deviation from the 
established theoretical predictions and has triggered the 
reconsideration of the theoretical description of these 
form factors (I10[7l[8l[9l[Il[IIl[I2[l3l[Il. Over the 
last year several theoretical advances appeared, which 
were able (or at least claimed to be able) to reconcile 
theory and experiment (e.g., (61Q3)) for the 7r°, 77 and K]' 
simultaneously. Nevertheless, there are also indications 
that the data of Ref . Q on n° cannot be easily accom- 
modated by theory (see Q [13 [HI III El EI). There " 
fore, there are still open questions and further research 
is foreseen both on the theory and the experiment side. 

The new experiments, e.g., KLOE-2 fT8l [T9lL and 
BES-III [ 20 1 will soon start measuring the two-photon 
form factors of the light pseudoscalars and related ob- 
servables. In case of a double-tag experiment, see 
Fig. [T] (right), the form factor as a function of both pho- 
ton virtualities is potentially of interest. Currently, the 
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double-tag experiments are the one performed with the 
KEDR detector at Novosibirsk [21 ] and the one per- 
formed by KLOE-2 Collaboration at Frascati lfT8l[T9l . 
Unfortunately, in the above experiments the shape of the 
form factor as a function of both photon virtualities is 
not accessible. Notice that the form factor unknown de- 
pendence on the second photon virtuality has possible 
impact on the precision of the experimental extraction 
of the quoted values of the form factor even in case of a 
single-tag experiment, as has been stressed in [22 ] and 
studied in l23l . In any case, a realistic simulation has to 
account for both photon virtualities. 

The transition form factor description is not only 
a milestone in our understanding of the U(l) axial 
anomaly [9], general principles of the strong interac- 
tion dynamics [11] and SU(3) flavor symmetry break- 
ing |24| . but it is also very important in the context of 
the Standard Model tests at low energy Q, mainly in 
the problem of the anomalous magnetic moment of the 
muon a^. The transition form factors (though in a dif- 
ferent kinematic region) are required for the calculation 
of the hadronic light-by-light (hLbyL) scattering part of 
the ^ usually denoted as ^ byL j^[^|27l[28ll291[30l. 
In view of the new g - 2 experiments to be carried out at 
Fermilab [ 31 ] and JPARC [32], it is crucial to improve 
the accuracy of a^ hyL ^ In context of hLbyL the im- 
provement of the models for the two-photon form fac- 
tors is important, which requires the high precision of 
the form factor measurements H9l . In turn, precision 
measurements require precision simulation. 



This subject has been a topic of the "INT Work- 
shop on the Hadronic Light-by-Light Contribution to the 
Muon Anomaly" (Seattle, February 28 - March 4, 2011) 
|http://wwwint.washingtoiT edu/PROGR AMS/ll-47w/| 
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Figure 1: The concept of single-tag (left) and double-tag (right) e + e —> e + e V experiment. 



Summing up, the state of the art is that both theory 
and experiment need a computer program for the pro- 
cess e + e~ — > e + e~ ( P- the Monte Carlo generator, which 
would allow to simulate the events under the condi- 
tions similar to the experimental ones. Our objective 
is to develop and to maintain a free and public tool for 
e + e~ — > e + e~ e P simulation (7r°, n, 77') to be used by ex- 
perimentalists and theorists. The guiding principles are 
the following: 

• to account for light meson phenomenology in de- 
tail via precise modeling of the form factors, in 
agreement with available data; 

• full description of both virtual photons; 

• high generator efficiency: accounting for the peak- 
ing behaviour of the cross section. 

This approach is being realized in the Monte Carlo gen- 
erator EKHARA (331 EH In the following Section we 
review the current status of this program. 

2. EKHARA Monte Carlo generator 

The Monte Carlo generator EKHARA is a computer 
program written in FORTRAN using the quadruple preci- 
sion numbers. The principal description of the program 
is given in the dedicated paper |33l . The generator can 
be used a stand-alone code or be interfaced to another 
(e.g., an experiment's) simulation program. 

The first version of EKHARA (33 El allows to sim- 
ulate the reaction e + e~ — > n + n~e + e~ with the focus on 
a specific event selection, namely that of the radiative 
return method for cross section cr(e + e~ —> mea- 
surement. It was used for the background (36l simu- 
lation by KLOE Collaboration 1571 1551 . This channel 
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can be important for the radiative return program of the 
BES-III experiment (39H20). The e + e~ -> 7r + 7Te + 6>- 
mode in EKHARA has the following features, which a 
user should be aware of: it generates weighted events 
and the applied phenomenology of the 7*7* — > n + n~ 
subprocess is not tuned to reflect the line shapes known 
from the dedicated 7*7* — > experiments. An im- 
provement of the above issue is on the agenda. 

The second version of EKHARA [33 ] features a new 
channel, e + e~ —> e + e~ ( P. The implementation of this 
channel and the applied phenomenology is oriented on 
the 7*7* physics. The phase space construction and the 
matrix element calculation is based on the exact formu- 
lae and the exact kinematics (in contrast to the Equiv- 
alent Photon Approximation, which was used in many 
generators in this field). 

It is important to stress that the (user-provided) form 
factors in EKHARA have to be meaningful both in the 
time-like and space-like regions of the photon, because 
both s- and ^-channel amplitudes and their interference 
can be switched on in the simulation, see Fig. [2] In gen- 
eral, the provided form factors have to be functions of 
two photon virtualities. The above criteria considerably 
reduces the user's choice, because the majority of the 
published formulae for the form factors within different 
theoretical approaches hold only for the case with one 
photon being real and the other — space-like and vir- 
tual. Strictly speaking, from the point of view of the 
numerical simulation, such a case is very specific and is 
never realized for the ^-channel amplitude because the 
leptons are not massless. 

In order to support a user, we have recently (23l sup- 
plied EKHARA with the transition form factors for 7r°, n 
and r{ mesons. Inclusion of r\ and n' meson produc- 
tion in EKHARA is a new feature with respect to (331 . 
The form factor modeling [23 ] was performed in the 
framework of the effective field theory, namely the chi- 
ral theory with resonances [40, 4l]|42) was used. The 
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The future plans of EKHARA development include the 
electromagnetic radiative corrections for the channel 
e + e~ — > e + e~ ( P. After the inclusion of the radiative cor- 
rections we plan to perform a detailed comparison of 
our generator with GGRESRC 1431 , which was used by 
the BaBar experiment (HE). The major update of the 
e + e~ — > e + e~n + 7r~ channel in context of possible y*y* 
physics motivation is also foreseen. 



Figure 2: The s- and f-channel amplitudes in e + e — > e + e P. 

quality of the model is good: by fitting only one free 
parameter to the data for 7r°, r] and r{ form factors (si- 
multaneously) we achieved x 2 /d.o.f. « 1.2. The exist- 
ing data from single-tag experiments in space-like re- 
gion (CELLO [3 ], CLEO [4] and BaBar ffl 0) were 
used not only as a reference for model improvement, 
but also for the generator verification. The verification 
is done by comparison of the simulated cross sections 
-^t for e + e~ — > e + e~P under the event selection simi- 
lar to that of the single-tag experiments. The pattern of 
data/simulation agreement is the same as for the form 
factors and this is an important consistency check. 

Recently, EKHARA was used by KLOE-2 Collabora- 
tion [19 ] for the feasibility studies of e + e~ — > e + e~ n°. 
On the basis of simulation, a per cent level of preci- 
sion is foreseen for the measurement of r^o^ rr by the 
KLOE-2 experiment at Frascati. KLOE-2 experiment 
can also perform the first measurement of pion transi- 
tion form factor F n o y * y (Q 2 ) in the space-like region in 
the vicinity of the origin: 0.01 < Q 2 < 0.1 GeV 2 . The 
impact of the proposed measurements on the evalua- 
tion of the Standard Model prediction for the anomalous 
magnetic moment of the muon, a^, was estimated lfT9l . 

3. Summary 

A brief review of the Monte Carlo event generator 
EKHARA is presented with the focus on the physics mo- 
tivation and several new features. It can be used for 
the simulation of the processes e + e~ — > e + e~ ( P and 
e + e~ — > e + e~7r + 7i~. The channels e + e~ — > e + e~7t°, 77, 
K]' are important for physics programs of the mod- 
ern experiments (e.g., KLOE-2 and BES III) and can 
be used for the simulation of the two-photon form fac- 
tor measurements at meson factories. The program is 
supplied with the form factors for 7r°, 77 and 77' mesons 
in agreement with the modern data on the transition 
form factors and corresponding differential cross sec- 
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